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Abstract 

We exhibit a specific implementation of the creation of geometrical phase through the 
state-space evolution generated by the dynamic quantum Zeno effect. That is, a system 
is guided through a closed loop in Hilbert space by means a sequence of closely spaced 
projections leading to a phase difference with respect to the original state. Our goal is the 
proposal of a specific experimental setup in which this phase could be created and observed. 
To this end we study the case of neutron spin, examine the practical aspects of realizing 
the "projections," and estimate the difference between the idealized projections and the 
experimental implementation. 



1 Introduction 



The effect of the observer in quantum mechanics is perhaps nowhere more dramatic 
than in the collection of phenomena loosely (and casually) known as the "quantum 
Zeno effect." This was first formulated by von Neumann P], |2|], and is deeply rooted 
in fundamental features of the temporal behavior of quantum systems ||. During 
the last decade there has been much interest in this issue, mainly because of an 
idea due to Cook [|J, who proposed using two- level systems to check this effect, and 
the subsequent experiment performed by Itano et al. [|5||. New experiments were 
proposed, based on the physics of the simplest of two-level systems: Neutron spin 
and photon polarization || 0. 

Most of the referenced papers deal with what might be called the "static" version 
of the quantum Zeno effect. However, the most striking action of the observer is not 
only to stop time evolution (e.g., by repeatedly checking if a system has decayed), 
but to guide it. In this article we will be concerned with a "dynamical" version of 
the phenomenon: we will show how guiding a system through a closed loop in its 
state space (projective Hilbert space) leads to a geometrical phase || || [10], [EL], |12| . 
This was predicted on general grounds 0], but here we use a specific implementation 
on a spin system |]I4[ and propose a particular experimental context in which to see 
this effect. It is remarkable that the Berry phase that is discussed below is due to 
measurements only: no Hamiltonian is needed. 

2 Forcing the pot to boil 

We summarize the main features of the quantum Zeno effect (QZE). Prepare a quan- 
tum system in some initial state ip(0). In time dt, by the Schrodinger equation, its 
phase changes by 0(dt) while the absolute value of its scalar product with the initial 
state changes by 0(dt 2 ). 

The dynamical quantum Zeno effect exploits the above features and forces the 
evolution in an arbitrary direction by a series of repeated measurements: Let ip 
evolve with the Hamiltonian H, so that in the absence of observations its evolution 
would be if)(T) = exp(—iHT)ip(0) (we take h = 1 throughout). Let there be a family 
of states 4>k, k = 0, 1, . . . , N, such that 0o = "^(0), and such that successive states 
differ little from one another (i.e., \{4>k+i\4>k)\ is nearly 1). Now let 5T = T/N and 
at Tf, = kST project the evolving wave function on 0^. Then for sufficiently large N, 
ip(T) « <p N . [The usual QZE is the special case 4>k = <po(— "0(0)) V k.] 

In the following we consider an experiment involving a neutron spin. It should 
be clear, however, that our proposal is valid for any system with the same two-level 
structure. 

2.1 Evolution with no Hamiltonian 

Assume first that there is no Hamiltonian acting on the system: one can think, for 
instance, of a neutron crossing a region where no magnetic field is present. The 
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time-evolution is due to measurement only. 

The system starts with spin up along the z-axis and is projected on the family of 
states 

/ 1 \ ak 
<f) k = exp(-i6 k (T ■ n) ( J with 9 k = — , k = 0, . . . , N , (2.1) 

where <r is the vector of the Pauli matrices and n = (n x , n y ,n z ) a unit vector (inde- 
pendent of A;). 

We assume that the system evolves for a time T with projections at times T k = k5T 
(k = 1, . . . ,N and 5T = T/N). The final state is 0o = (o) 

|V(T)) = \<Pn)(^N-i) ■ ■ ■ (0 2 |0l)(0l|0 O ) 
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exp(mn z )|0jv) 

= exp(mn 2 ) exp(— mtr • n)|0 o )- (2-2) 

Therefore, as — > oo, ip(T) is an eigenfunction of the final projection operator Pn, 
with unit norm. If cos = n z and a = n, 

ip(T) = exp(i7r cos 0)(— l)0o = exp[— m(l — cos9)]0o = exp(— iQ/2)(p , (2.3) 

where Q is the solid angle subtended by the curve traced by the spin during its 
evolution. The factor exp(— ifl/2) is a Berry phase and it is due only to measurements 
(the Hamiltonian is zero). Notice that no Berry phase appears in the usual quantum 
Zeno context, namely when <p k oc <po V k, because in that case a = in (|2.2[). 



To provide experimental implementation of the mathematical process just de- 
scribed, one could (in principle) let a neutron spin evolve in a field-free region of 
space. With no further tinkering, the spin state would not change. However, suppose 
we place spin filters sequentially projecting the neutron spin onto the states of Eq. 
( |2.1| ), for k = 0, . . . , N. Thus the neutron spin is forced to follow another trajectory 
in spin space. The essence of the mathematical demonstration just provided is that 
while N measurements are performed, the norm of wave function that is absorbed by 
the filters is iV-0(l/iV 2 ) =0(1/ N). For iV — >• oo, this loss is negligible. Meanwhile, 
as a result of these projections, the trajectory of the spin (in its space) is a cone 
whose symmetry axis is n. By suitably matching the parameters, the spin state can 
be forced back to its initial state after time T (T3j . 



It is interesting to look at the process ( |2.2j ) for iV finite. The spin goes back to its 
initial state after describing a regular polygon on the Poincare sphere, as in Figure 
la. After A^(< oo) projections the final state is 



a . . a- N 



\ip(T)} = ( cos — + m z sin— ) exp(-iacr • n)|0 o ). (2.4) 
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a 

Figure 1: a) Spin evolution due to N = 



b 

5 measurements, b) Solid angles. 



For a = it the spin describes a closed path and 
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\ip(T)) = (cos— + m 2 sin — ) exp(-ivr)|0 o ) 



N_ 

7T r, 7T\2 / . „ / vr 



cos — + n z sin —J exp \iN arctan f n^tan — J J exp(— Z7r) | 0o) - 

(2.5) 

The first factor in the far r.h.s. accounts for the probability loss (N is finite and there 
is no QZE). We can rewrite Q2.5Q in the following form 

|^(T))= PiV exp(-^)|0o), (2.6) 

where 

Pn = [cos — + n z sm—\ , (2.7) 
(3n = 7T — N arctan ^cos tan . (2.8) 
In the "continuous measurement" limit (QZE), we have 

P = J im PN = 1, 

/3 = lim /3iv = 7r(l-cose) = ^, (2.9) 
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where Q is the solid angle subtended by the circular path, viewed at an angle 6 (see 
Figure la). We recover therefore the result (p!3|). 

The relation between the solid angle and the geometrical phase is valid also with 
a finite number of polarizers N. Indeed, it is straightforward to show that the solid 
angle subtended by an isosceles triangle with vertex angle equal to 2a (Figure lb) 
has the value 

VL 2a = 2a — 2 arctan(cos tana). (2-10) 

Hence if the polarizers are equally rotated of an angle 27r/N, the spin describes a 
regular iV-sided polygon, whose solid angle is 

tt {N) = Ntt 27r/N = 2rr-2N arctan (cos 9 tan ^\ = 2f3 N , (2.11) 

where we used the definition ( |2.8|) . This result is of course in agreement with other 
analyses based on the Pancharatnam connection ||. 

The above conclusion can be further generalized to the general case of an arbitrary 
(not necessarily regular) polygon. Indeed, if the polarizers are rotated at (relative) 
angles a n with n = 0, . . . , N, so that 

JV 

E 2 «n = 27r ' ( 2 - 12 ) 
n=l 

the solid angle is 

N N 

f2/jy\ = ^2 ^2a„ = 27T — 2 arctan(cos0tana n ). (2-13) 



n=l n=l 

This is also twice the Berry phase. Notice that if all a n — > as N — ► oo one again 
obtains the limit ( |2.3| ): 

N 

Q' = lim Q' N = 2tt - 2 lim V a n cos 9 = 0. (2.14) 

n=l 

We emphasize that these predictions for the N < oo case are not trivial from the 
physical point of view. The above phases are computed by assuming that, during a 
"projection" a la von Neumann, the spin follows a geodesies on the Poincare sphere. 
The mathematics of the projection has no such assumptions. The "postulate's" only 
job is to relate all this projection formalism to measurements. 

2.2 Evolution with a non-zero Hamiltonian 

Let us now consider the effect of a non-zero Hamiltonian 

H = fiab, (2.15) 

where b = (b x , b y , b z ) is a unit vector, in general different from n. One can think of a 
neutron spin in a magnetic field. See Figure 2. 
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If the system starts with spin up it would have the following — undisturbed 
evolution: 



ip(t) = exp(— ifxta ■ b) 



(2.16) 



Now let the system evolve for a time T with projections at times T k = k5T (k 
1,...,N and 5T = T/N) and Hamiltonian evolution in between. Defining Pq 



1 00 ) (00 
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0/ 



the 2x2 projection operator at stage-fc is 

Pk = \4>k){4>k\ = exp(-i8 k cr ■ n)P exp(i8 k cr ■ n) 
and the state evolves to 



V>(T) 



N 



J] [P k exp(-in5Ttr-b)] 



Lk=l 



(2.17) 



(2.18) 



where here and in subsequent expressions a time-ordered product is understood [with 
earlier times (lower k) to the right]. Using Pq = Pq, Eq. ( |2.18| ) can be rewritten 



ip(T) = exp(— iacr ■ n) 



N 



n bh 

.k=l 



with 



B k = Pq exp(z#fc£r ■ n) exp(—in5Tcr ■ b) exp(—i9h-i<r ■ ti)Pq 



(2.19) 



(2.20) 
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(9q = 0). The computation of B k requires a bit of SU(2) manipulation. By using 

[a A,(T B] = 2ia ■ Ax B (2.21) 
(<r ■ A)(cr • B)(cr ■ A) = 2(A ■ B)cr ■ A — (A- A)cr ■ B, (2.22) 

valid for c-number A and B, one gets 

exp(z$cr ■ n)er ■ bexp(—ider ■ n) — cr ■ b, (2.23) 

with 

6(0) = b cos 29 + n(b ■ n) (1 - cos 26) + b x n sin 26, (2.24) 

which is the vector b rotated by 26 about the n-axis. The calculation of B k is now 
straightforward: 

Bk = Po exp(i56er ■ n) exp(—ifi5Ter ■ b(6k-i))Po 

= P (l+i66(T -n-itiSTtr -b(6 k )) P + O(l/N 2 ), (2.25) 

where 56 = 6^. +1 — 6^ is fc-independent. Second order terms in 1/N drop out when 
the product fl2.19p is computed for iV — > oo, so that 



N 

n B k 

k=l 



N 

[] Po(l + iS6a ■ n - ifxSTa ■ b(6 k ))P 

k=l 
N 

]J {^o + iP (58<r ■ n - fiSTcr ■ b(6 k ))P } 

k=l 
N 

n^o{l + i[^*-^6^ fc )]} 

k=l 

{ N 1 
P expliY^ (&0n z - jj5Tb z (6 k ) 



k=l 



(2.26) 



where we have used P a x P = P a y P = and P <j z P = P . The continuum limit 
can be computed by letting the summations in ( |2.26 ) become integrals in dT and dd. 
Moreover, ^ = — , which enables one to change integration variable and get for the 
"(1,1)" component of Y[ k= i B k (all other components being zero) 

exp (in z [ dd — ifi— [ [b z cos 26 + (6 • n)n z {l — cos 26) + (6 x n) z sin 20] dd 
\ Jo a Jo 



T 

exp [ in z a — in— 
a 



sin 2a 



+ (b ■ n)n z a - 



sin 2a 



+ (b x n) 



1 — cos 2a 



(2.27) 



The final state is an eigenstate of Pjy with unit norm, independent of the Hamiltonian 
H: 



T 

exp I —in— 

a 



sin 2a 



x exp [%an z — iacr ■ n) 



+ (b ■ n)n z ( a — 

1 




sin 2a 



f b x n) 



1 — cos 2a 



(2.28) 
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The first factor in ( 2.28|) is obviously the "dynamical phase." Note that up to a phase, 
if)(t) is just 0^, with k = tN/T. Therefore 



T r a 

(ip(t)\H\ip(t))dt = — ((fi \exp(i8cr ■ n)ficr ■ bex-p(—i6cr ■ n)\(fi )dO 

a Jo 



sin 2a ( sin 2a 

Or.— h o • n)n z 1 - 



2a 



2a 



[b x n) 



1 — cos 2a 



2a 



(2.29) 



because the phases drop out in the above sandwich. It follows that the remaining 
phase in ( p.28|) , when the spin goes back to its initial state, is the geometrical phase. 
When a = 7i 

4>(T) = exp (—ifi/2) exp (-i/xT(6 • n)n z ) \\ , (2.30) 

where fl is the solid angle subtended by the curve traced out by the spin, as in ( j2.3|) , 
and fj,T(b-n)n z yields the dynamical phase, as can also be seen by direct computation 
of ( |2.29|) . We remark that if time ordered products are looked upon as path integrals 
|fL6[| ) then our above demonstration is effectively a path integral derivation of the 
geometrical phase. 

A practical implementation of the process just described would involve an exper- 
imental setup similar to the one described after Eq. (|2.3| ), but with a magnetic field 
whose action on the spin is described by the Hamiltonian (|2.15 ). If the neutron were 
to evolve only under the action of the Hamiltonian, its spin would precess around the 
magnetic field. However, the sequence of spin filters, which project the neutron spin 
onto the states (|2.1|), compel the spin to follow the same trajectory as in the previous 
case [Eq. ( |2.2j )], i.e. a cone whose symmetry axis is n. As above, the spin acquires a 
geometrical phase, but now there is a dynamical phase as well. 



2.3 A particular case 



It is instructive to look at a particular case of ( |2.28|) - (|2.30|) . We first note that if 



/i = in (|2.28|) we recover (|2.2|) . Now let b = n. In this situation the projectors and 



the Hamiltonian yield the same trajectory in spin space (although, as will be seen, 
at different rates). If /x = (so that H = 0), the spin evolution is only due to the 
projectors and the final result was computed in 



ip(T) = exp(-ifi/2)0 O - (2.31) 



If, on the other hand, there is a nonvanishing Hamiltonian (|2.15|) , but no projectors 
are present, a cyclic evolution of the spin is obtained for \iT = ix. The calculation is 
elementary and yields 

ip(T) = exp(-m)0 o - (2.32) 
Observe that the dynamical phase in this case is [/iT = n,b = n and a = n in Eq. 

(Pi] 



T 

(ip(t)\H\ip(t))dt = 7m z = tt[1 - (1 - n*)] = 7T - Q/2. (2.33) 
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Therefore, the 'V phase in ( J2.32[ ) can be viewed, a la Aharonov and Anandan |I3~ 
as the sum of a geometrical (0/2) and a dynamical (n — Q/2) contribution. 



Now let both the Hamiltonian and the projectors be present. From Eq. Q2.30| ). 

one gets 

i/>(T) = exp {-in/2) exp (-ifiTn z ) Q , (2.34) 

Notice that the value of \i is now arbitrary, so that \iT is not necessarily equal to ix 
(the cyclic evolution of the spin is due to the projectors, not to the Hamiltonian). 
When /xT < ir, the projections are too "fast" and do not yield ( |2.32|) . On the other 
hand, when \xT > it, the projections are too slow and supply less phase, in comparison 
with Eq. (|2.32|) . Only in the case pT = it do the projections yield the right phase 
in ( |2.32j ). Their presence is superfluous in this case: one would obtain exactly the 
same vector and the same phase without them. Our conclusions are summarized 
in Table 1. In some sense, one may say that the Hamiltonian dynamics provides a 
"natural clock" for the phase of the wave function. 

Table 1: Phases for cyclic spin evolutions 





H = 


H = jicr ■ b 


H = jicr ■ b 




and projections 


no projections 


and projections 


0geom 


0/2 


0/2 


0/2 


0dyn 





7i - 0/2 


fiTn z 


0tot = 0geom + 0dyn 


0/2 


7T(= flT) 


0/2 + pTn z 




cyclic evolution 


cyclic evolution 


cyclic evolution 




due to projections 


due to field 


due to projections 



3 A Gedanken Experiment 

An experimental implementation with neutrons would be difficult because it would 
involve putting a QZE set-up inside an interferometer in order to measure phase. We 
therefore restrict ourselves to a gedanken experiment based on the use of 3 He as a 
neutron polarization filter ||17|| . It is well known Jl8| that Helium 3 
neutrons but polarized 3 He only absorbs one spin state of a neutron beam- 



'black" to 
-hence 

acts as a 50% absorber of a beam; the rest of it emerges fully polarized. In practice 
an external magnetic field is used to maintain the polarization axis of the 3 He. If this 
external bias field were to be given a slow twist along a longitudinal axis, the state 
of polarization of the 3 He should follow the direction of the twist. A neutron beam 
propagating through a cell of high-pressure polarized 3 He along an axis aligned with 
the direction of twist will become fully polarized and should develop a Berry phase 
according to the argument of the previous section. 

From an experimental perspective a significant problem is that we so far lack a 
notion of slowness (as when we speak of "slow twist" of the B field). In the previous 
calculation, it is implicitly assumed that 9 changes more slowly than t (time): in other 
words, the relaxation processes in the 3 He are given enough time (are fast enough) to 
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function as a polarizer. A full treatment of this problem should therefore describe the 
physics of the projection process. We now tackle this issue and see that the notion of 
slowness can be given quantitative meaning in terms of a condition for adiabaticity. 

In practice, the absorption of the non-selected spin state occurs over a finite 
distance, of the order of one or two centimeters. This situation can be modeled 
via the following family of effective (nonhermitian) Hamiltonians: 



H k = -iV\<t>i)(<t>t\, 



where V is a real constant and 
<pj: = exp(—i9 k cr ■ n) 



with 9 k = — , k = 0, 



N . 



(3.1) 



(3.2) 



Note that {4> k \4> k ) = [ see Eq. flO)]. We first assume, for simplicity, that no external 
( 3 He aligning) magnetic field is present. We define 

Ph = I0*><0* I = exp(-Mk<r • n)J^ exp(tf*ir ■ n) (P L = . (3.3) 

Obviously P k = 1 — P k , where P k was defined in (|2.17|) . The evolution engendered 
by the above Hamiltonian reads 



-iHuT 



P k + eP k = exp(-i6 k cr ■ n) 



exp(i9 k a ■ n) = P' k , 



where (inserting h) 



e = e- VT / h 



(3.4) 



(3.5) 



is a parameter yielding an estimate of the efficiency of the polarizer. One can estimate 
a minimal value for V: for a thermal neutron (speed v ~ 2000m/s) and an absorption 
length i on the order of 1 cm for the wrong-spin component, one gets r = tjv ~ 5/xs 
and one obtains a good polarizer for V > H/r ~ 10~ 29 J ~ 10~ 7 meV. 
The evolution can be computed by using the technique of Section 

Po + ^p^y. 

- N 



ip'(T) = exp(—iacr ■ n) 



with T = Nt and 

N 

I[B> = 

k=l 



.k=l 



N 



(3.6) 



N 

II s/nil + iSea-n^P^ n p o 

k=l 

1 + i59n z i89e 1/2 n_ 
efl — i59n. 



iJp^(S9cr ■ n) 



k=l 



me 1 / 2 



N 



where n± = n x ± in y . The evaluation of the above matrix product when N 
lengthy but straightforward. One gets 

ip'(T) = exp(—iaa- ■ n)A4(p , 



(3.7) 



OO IS 



(3i 
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where 



_ / A ch(aA) + (b + in*) sh(aA) in_ sh(aA) \ , . 

A V zn+ sh(aA) A ch(aA) - {b + in*) sh(aA) J ' ^ ^ 

with 

. VT 



A = y& 2 + 2ibn z - 1. (3.10) 

We are interested in the limit of large b = VT j2ah. Indeed, larger values of b 
correspond to more ideal polarizers. In fact 7 = V/h represents the absorption rate 
of the wrong component of the spin, while uo = 2a /T is the angular velocity of 
precession (the spin describes an angle of 2a in time T). The parameter b = 7/0; is 
the ratio of these two quantities. Large values of b imply 

7>^, (3.11) 

i.e., an absorption rate much larger than the velocity of precession. In other words, the 
spin rotation must be sufficiently slow to allow the absorption of the wrong component 
of the spin. By introducing the neutron speed v , one can define the absorption length 
I = u/7 = vh/V and the length covered by the neutron while rotating for 1 rad, 
L = v/uj = vT/2a. Hence (|3.11| ) reads 



L > L (3.12) 



These are all conditions of adiabaticity. 

In the large b limit, using the definition Q3.10|) , ( |3.9|) becomes 

e a(A-b) / A + & + 



1 1— n? -n- \ /I 

1 — a—^ Zttt \ , ^ / 1 




2b 2b 

; 1 - w 



O - . (3.13) 



Remembering the definition of b in ( |3.10| ), one gets 

— > exp(ian z )Po, when ^— — > 00. (3.14) 

The above formula yields the first corrections to an ideal, purely adiabatic evolution. 
Basically, the system is projected on slightly different directions, thereby rotating in 
spin space. But if the system "on its own" (i.e., through its dynamics) manages to 
rotate significantly between projections, then more will be absorbed on the next pro- 
jection and it will not follow the rotating field, at least not without loss of probability 
(or intensity). 
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It is interesting to note that the same result can be obtained by considering a 
continuous version of the effective Hamiltonian (13.11) 



H{t) = -iVP L {t) = -iVU ] {t)P^U{t), (3.15) 

where 

U(t) = exp (i^t a ■ nj (3.16) 

is a unitary operator (rotation). The state vector ip(t) satisfies the Schrodinger equa- 
tion 

idt4>(t) = H(t)1>(t). (3.17) 
Consider now the following rotated vector 

4>{t) = U(t)m- (3-18) 

It is easy to prove that it satisfies the equation 

idt$(t) = H$(t), (3.19) 

where 

H = iU{t)U ] {t) + U{t)H{t)U ] {t) = ~tr ■ n - iVP^ (3.20) 
is independent of t. One then gets 

ijj{t) = U\t)4>{t) = exp (-i^t ° ■ nj exp(-iHt)ip(0), (3.21) 

where 



HT=-a,.n-iVTP- = -aM, M = ^ _ n ^ +i2b ), (3-22) 

b being defined in (|3.10p . Hence one obtains 

exp{-iHT) = exp(iaM) = M (3.23) 
and ( gg| ) yields (p|). Observe that 

H = -u ijPj~, (3.24) 



from which it is apparent the previous interpretation of the coefficients u and 7. 

The above calculation was performed by assuming that no external field is present. 
However, we do need an external B field, in order to align 3 He. Its effect can be 
readily taken into account by noticing that, when the neutron crosses the region 
containing polarized 3 He, if the conditions for adiabaticity are satisfied, the neutron 
spin will always be (almost) parallel to the direction of 3 He and therefore to the 
direction of the magnetic field. The resulting dynamical phase is therefore trivial 
to compute and reads 0d yn — nBT/K. In order to obtain the geometric phase in 
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a realistic experiment, such a dynamical phase should be subtracted from the total 
phase acquired by the neutron during its interaction with 3 He. Incidentally, notice 
that this is experimentally feasible: one can take into account the contribution of a 
large dynamical phase due to the magnetic field and neatly extract a small Berry 
phase [19|. The novelty of our proposal consists in the introduction of polarizing 3 He 
to force the neutron spin to follow a given trajectory is spin space. 

An alternative realization relies on a set of discrete 3 He polarization filters with 
progressively tilted polarization axes, as a finite-difference approximation to the sys- 
tem discussed above. Such a system would be a neutron analog of a set of polaroid 
filters with progressively tilted axes through which a photon beam propagates with 
little or no loss (in the limit of small angles) as proposed by Peres [0. However, in 
the case discussed in this Letter, the axes of the neutron polarizers need not belong 
to a single plane and the neutron can acquire a Berry phase as well as change in 
polarization direction. 
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